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ABSTRACT 

Testing Newtonian gravity in the weak-acceleration regime is vital to our understand- 
ing of the nature of the gravitational interaction. It has recently been claimed that the 
velocity dispersion profiles of several globular clusters flatten out at large radii, remi- 
niscent of galaxy rotation curves, even though globular clusters are thought to contain 
little or no dark matter. We investigate this claim, using AAOmega observations of 
four globular clusters, namely M22, M30, M53 and M68. M30, one such cluster that 
has had this claim made for its velocity dispersion, was included for comparison with 
previous studies. We find no statistically significant flattening of the velocity disper- 
sion at large radii for any of our target clusters and therefore we infer the observed 
dynamics do not require that globular clusters are dark matter dominated, or a modi- 
fication of gravity. Furthermore, by applying a simple dynamical model we determine 
the radial mass-to-light profiles for each cluster. The isothermal rotations of each clus- 
ter are also measured, with M22 exhibiting clear rotation, M68 possible rotation and 
M30 and M53 lacking any rotation, within the uncertainties. 
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1 INTRODUCTION 

Newtonian gravitation has been shown to describe the mo- 
tions of bodies with intermediate accelerations (e.g. Solar 
System bodies) very accurately. The Newtonian theory of 
gravity breaks down in the presence of strong gravitational 
fields, where it has been successfully superseded by Gen- 
eral Relativity, but there is evidence that there are also dis- 
crepancies in the low-acceleration regime. Spiral galaxies, 
for example, require the invocation of Dark Matter (DM) 
to reconcile the discrepancy between their rotational ve- 
locities and our understanding of Newtonian gravity. An- 
other indication is the so-called "Pioneer anomaly": radio- 
metric data from Pioneers 10 and 11 have shown that both 
are experiencing an unexplained constant accele ration of 
a = ( 8.74±1.33) x 1CT 10 ms~ 2 toward the Sun (see lde Diegol 
120081 . and references therein). In fact, all spacecraft in the 
outer Solar Syst em seem to be experien cing this anomalous 
acceleration (e.g. lAnderson et al. 12 002). potentially indicat- 
ing a deviation from Newtonian gravity at low accelerations. 
It has been claimed that one of the leading versions of 
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modified Newtonian dynamics (MOND; Milgrom 1983|) can 
predict the kinematic properties of galaxies without invok- 
ing DM. MOND predicts a breakdown of Newtonian gr avity 
at an acceleration scale a » 1.2 x 10~ 10 ms" 2 (see lloriol 



120091 . and references therein). This is approximately the ac- 
celeration regime where DM becomes necessary to reconcile 
theoretical velocity profiles of galaxies with observation. 

Globular clusters (GCs) are the ideal test bed for 
this divergence from Newtonian predictions, since they are 
thought to contain little or no D M. This is base d on evi- 
dence from dyna mical models (e.g. |Phinnevlll993h . -/V-body 
simulation s fe.g. iMoorel Il996l), ob servations of GC tidal 



tails (e.g. lOdenkirchen et al. I 200l[). dynamical a nd lumi- 
nous masses of GCs fe.g. iMandushev et al. |[l99ll ) and the 
lack of microlensing events from GC-ma ss dark haloes (e.g. 
iNavarro rf"anil997l ; llbata et al. Il2002t ). In addition, GCs 
at different distances from the Galactic centre experience 
differing gravitational attractions from the Galaxy, so that 
any Galactic influence can be ruled out if all exhibit sim- 
ilar behaviour. Furthermore, the accelerations experienced 
by s tars in GCs drop b elow ao well inside the tidal radius 
fe.g. lScarpa et aZ.ll2007l). 

Recently, Scarpa et al\ l|2007l ) showed that there may 
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be a flattening of the velocity profiles of several GCs in the 
Galactic halo (w Centauri, M15, M30, M92 and M107), at 
or about the acceleration where DM is invoked to main- 
tain stability i n rot ating galaxie s. On the other h and, 
iMoffat fc Tothl l|2008h overlaid the IScarpa et al\ (|2007f ) re- 
sults with a Modified Gravity (MOG) model and found lit- 
tle, or no, deviation from Newtonian gravity, for GCs with 
masses less than a few times 1O 6 M0. Because both groups 
used an identical set of measured velocity dispersions, there 
is an obvious discrepancy. If MOND were correct, it would 
mean that all gravitational interactions would diverge from 
Newtonian gravitation in the predicted regime. Irrespective 
of MOND itself, testing the gravitational interaction at low 
accelerations is extremely important in the overall under- 
standing of gravity. This paper continues this test, with four 
Galactic halo GCs, namely M22 (NGC 6656), M30 (NGC 
7099), M53 (NGC 5024) and M68 (NGC 4590). 



2 DATA ACQUISITION AND REDUCTION 

We used the multi-object, double-beam spectrograph 
(AAOmega) on the 3.9 m Anglo- Australian Telescope at Sid- 
ing Spring Observatory in New South Wales, Australia, to 
obtain the spectra for this survey. AAOmega has 392 fibres 
covering a 2 square degree field of view with each fibre capa- 
ble of obtaining a single spectrum from a single star in this 
field. The 1500V grating (with a resolution of R = 3700) was 
used in the blue arm and the 1700D grating (R = 10000) 
was used in the red arm, with the central wavelengths set to 
5200A and 8640A, respectively. This configuration was cho- 
sen to include the calcium triplet lines at 8498A, 8542A and 
8662 A, for accurate velocity determination and proxy metal- 
licities, as well as the swathe of iron and magnesium lines 
around 5200A for accurate metallicity measurements. The 
astrometric positions of the targets, taken from the 2MASS 
Point Source Catalogue^, have accuracies of ~ 0.1". 

The observations were taken over two observing runs, on 
5-8 and 14-19 June 2008, with an average seeing of ~ 1.5". 
We obtained 3x 1200-second exposures per fibre configura- 
tion, to obtain S/N~ 20 — 100, with several configurations 
per cluster (see Table [1] for actual numbers of spectra ob- 
tained for each cluster - note that for M68, 3 configurations 
were obtained in June 2008, and 5 further configurations 
were obtained during a separate service run in May 2009). 
In this study we also made use of data for M30 obtained 
in a previous ob serving run in 2006 (originally published by 
iKiss et aLl 120071 ). Flat fields, arc lamp exposures and ~25 
dedicated sky fibres were used for data reduction and cal- 
ibration. FeAr, CuAr, CuHe, CuNe and ThAr arc lamps 
were used to ensure accurate wavelength calibration. Data 
reduction was performed using the 2df dr pipelinfl designed 
specifically for the reduction of AAOmega data. The efficacy 
of the pipeline has been checked with a comparison of indi- 
vidual stellar spectra. 



Table 1. Clusters selected for the current survey, based on the 
criteria outlined in the text, with the total number of fibre con- 
figurations, and spectra, obtained during the run. Note that the 
number of spectra quoted here is the total number of spectra ob- 
tained before Galactic contaminants were removed. For M68, 3 
configurations were observed on the two observing runs in June 
2008 and 5 subsequent configurations were observed during a ser- 
vice run in May 2009. 



Cluster 


# Fibre Configurations 


Total # Spectra 


M22 


10 


3407 


M30 


2 


620 


M53 


6 


1727 


M68 


3+5 


2650 



2.1 Cluster Selection 

We aimed to observe several hundred stars from each of 
four GCs, namely M22 (NGC 6656), M30 (NGC 7099), M53 
(NGC 5024) and M68 (NGC 4590). These four targets were 
selected by several criteria: they are bright (My < —7), 
nearby (D < 30kpc) and have radial velocities > 1.5 a from 
the peak velocity given by the Besangon Galaxy Modejf] in 
that direction, for ease of extracting the cluster members 
from Galactic contaminants. One cluster (M30) was included 
specifi cally because it was one of the targets of lScarpa et all 
and it satisfied all other selection criteria. This allows 
a comparison with the results of that study. 

We produced a Colour Magnitude Diagram (CMD) for 
each cluster using data from the 2MASS Point Source Cat- 
alogue, and the Red Giant Branch (RGB) was identified in 
colour-magnitude space. Selection of individual target stars 
was based on the J — K colours and K magnitudes of the 
RGB of the cluster to minimise the number of stars selected 
from the Galactic population. Despite this selection process, 
a number of Galactic contaminant stars were still expected 
in the final sample and it was therefore necessary to remove 
them before analysis. Section [2.1.11 describes this process. 



2.1.1 Cluster Membership 

To select cluster members, it was necessary to determine sev- 
eral quantities for each observed star. Atmospheric parame- 
ters and radial velocities were determined using an iterative 
pro cess, combining best fits to the synthetic spectra from 
the iMunari et all ( 2005h spectrum library, with \ 2 fitting, 
and cross-correlating this best fit model with the observed 
spectrum to calculate the radial velocity. This approach is 
very similar to that ado pted by the Radial Velocity Experi- 
ment (RAVE) project (|Steinmetz et qZ.|[2006l ; IZwitter et all 
2008), an d is based on the same synthetic spectral library 
as RAVE. lKiss et all l|2007h outlined this method in detail. 

For M22, M30 and M68, we used four parameters to 
identify members, namely metallicity ([m/H]), equivalent 
width of the calcium triplet lines, radial velocity and the dis- 
tance from the cent re of the clust er [terminating at the tidal 
radius as quoted by Irlarrisl (|l996l )]; only stars that matched 
using all criteria were judged to be members. Figure Q] shows 
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all the stars for which spectra were obtained, highlighting 
those determined to be cluster members. Note that spectra 
from observing runs over 2006-2008 were combined with the 
spectra of M22, M30 and M68 taken during the current sur- 
vey to ensure greater statistical significance of the results. 

M53, being closer to the Galaxy in both [m/H] and 
radial velocity, required a non-degenerate metallicity to 
determine memb ership. We used the method outlined by 
ICole et ~al\ l|2004h to calculate [Fe/H] from the equivalent 
widths of the calcium triplet lines for each star, resulting in 
a much cleaner selection for this cluster. Further cuts were 
made on temperature (T e a < 9000 K) and the quality of the 
synthetic spectral fit to produce the final sample used for 
analysis. Typical uncertainties in the radial velocities of our 
member stars are ~ 3kms _1 , ~ 2kms _1 , ~ 3kms _1 and 
~ 2.5 kms -1 , for M22, M30, M53 and M68, respectively. All 
stars have uncertainties less than 6kms . 

The sample consisted of 345, 194, 180 and 123 RGB 
member stars for M22, M30, M53 and M68 respectively. 
Note that for M30, the 168 Main Seque nce Turn Off (MSTO) 
and sub-giant branch stars from the IScarpa et all (|2007l ) 
study were also included in our analysis, bringing the total 
number of members for that cluster to 362. The consequence 
of mixing stellar types in our analysis of M30 is discussed in 
Section EH1 



3 RESULTS 
3.1 Rotation 

The radial velocities were used to determine the rotation of 
each cluster, assuming an isothermal rotation. The rotation 
was measured by dividing the cluster in half at a given posi- 
tion angle (PA) and calculating difference between the aver- 
age velocities in the two halves. This was performed over PAs 
in steps of 10° for all clusters, except M53 was done in 15° 
steps to avoid aliasing effects. The resulting data were then 
fit with a sine function. Figure [2] shows the measured rota- 
tion for each c luster and the best-fit sine function. Note that 
IScarpa et all HmzI) detected no rotation in M30 to a level 
of 0.75kms _1 and we corroborate their result; we detect no 
evidence for rotation to a level of 0.8 km s -1 . For M53, no 
rotation is evident to 0.5 kms -1 , however, it is clear that 
M22 is rotating, with a projected axis of rotation approxi- 
mately North-South. M68 shows some indication of rotation, 
although additional data are required to confirm this. 



3.2 Velocity Dispersions 

To determine the velocity dispersions of our samples, we first 
corrected all four clusters for rotation, and then binned the 
stars by distance from the cluster centre (annuli), ensuring 
approximately equal numbers of stars per bin. The mean 
and standard deviation were calculated for the velocities in 
each bin. Our velocity dispersions were determined using 
a M arkov Chain M onte Carlo maximum likelihood method 
fsee lGregorvll2005l for an overview of MCMC methodology) , 
which takes into account the individual velocity uncertain- 
ties on the stars and provides the velocity dispersion in each 
bin with associated uncertainties. 



The resulting velocity disp ersion profi l es wer e then over- 
plotted with the best-fitting IPlumme r] (|l91ll ) model, a.s 
shown in Figure [3] The data points have been plotted at the 
mean radius of the stars in each bin. Because fewer cluster 
members exist at larger radii, the outermost bin is generally 
much larger than the inner bins. To ensure our Plummer 
model fits do not exhibit an artificial slope in the outskirts 
of the clusters, masking any possible flattening of the veloc- 
ity dispersion profiles, we also performed fits to the profiles 
with the last bin excluded. The fits are nearly identical for 
M22, M30 and M68. For M53 a small discrepancy between 
the fits was observed due to the scatter in the velocity dis- 
persions of the inner bins, however, this discrepancy has no 
impact on our overall results. The reduced \ 2 values for the 
Plummer fits shown in Figure|3]are ~ 2.3, ~ 1.1, ~ 2.0 and 
~ 1.5 for M22, M30, M53 and M68, respectively. We used 
two parameters in the fitting process, namely the central ve- 
locity dispersion and the scale radius (r s ; containing half the 
cluster mass). Due to the nature of the model it is possible 
to calculate the total cluster mass from the central veloc- 
ity dispersion (ctq) and r 3 fsee lDe ionghc 1987j for a detailed 
discussion of Plummer models and their application) : 



The total masses are shown in each panel of Figure [3] 

We chose the Plummer model to fit to our velocity 
dispersion profiles because it is a monotonically decreas- 
ing function and, therefore, any flattening of the profiles 
would be clearly s hown. The cluster in common between this 
study and that of lScarpa et all (|2007T ). M30, shows the best 
Plummer fit of the four. Note th at our complete da taset, 
which includes those data used bv lScarpa et all (|2007T l. has 
been binned to ensure similar numbers of star s per bin. Be- 
cause our binning is different to that used by IScarpa et all 
(2007), a data point-by-data point comparison between the 
two studies is not possible here. Within the limits of the 
model, and the uncertainties in the data, the velocity dis- 
persion does not exhibit any flattening at large radii, and 
therefore there is no indication that either DM or a mod- 
ified theory of gravity is required to explain the velocity 
dispersion of M30. The same case can be made for the other 
three clusters we examined, and the velocity dispersion pro- 
files, along with best fit Plummer models (and parameters) 
are presented in Figure [3] 

Different stellar populations in GCs are known to have 
differing velocity scatters with i n the cluster ("velo city jit- 
ter" ; see iGunn fc Griffin! H979I ; ICarnev et al1l2003l . for de- 
tailed discussions). This jitter is due to instabilities in the 
atmospheres of luminous giant stars. S ince we included the 
MSTO and sub-giant stars from the IScarpa et all (|2007l ) 
study with our RGB sample, this jitter should be evidenced 
as a slightly increased velocity dispersion from the result 
by IScarpa etail (|2007f l. This increase in dispersion will be 
smaller, however, than if we produced our velocity disper- 
sions from our RGB stars alon e. Comparing our d ispersion 
profiles for M30 with those of lScarpa et all l|2007l ). we see 
some evidence for this increase in the inner regions of the 
cluster, although the difference is well within the uncertain- 
ties. No incr e ase in dispersion is seen in the outer parts; 
IScarpa et all (|2007l ) quote a dispersion of ~2.25kms 1 at 
a radius of 25 pc, whereas our data show a dispersion of 
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Figure 1. Members of each cluster used for analysis (circled dots), based on the selection method described in the text. The uncircled 
dots are the stars which were observed and determined not to be cluster members. The large circle is the tidal radius of the cluster from 
lHarrisl(|l996h . In each panel, North is up and East is to the left. 
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Position Angle (N = 0. E = 90) Position Angle (N = 0, E = 90) 

Figure 2. The rotation of each cluster calculated as the difference between the mean velocities from each side of the cluster along equal 
position angles, as described in the text. The best x 2 fit sme function is overplotted, and a typical error bar is represented in the lower 
right of each panel. 



~1.6kms 1 at this distance, although this discrepancy is 
still within the uncertainties. 



The main result from the lScarpa et al. (2007) study was 
an apparent flattening of the velocity dispersion profiles of 
five of their six clusters, indicating a DM component, or 
a nee d for a modified grav ity theory to explain their re- 
su lts. iMoffat fc Tothl (l2008h overlaid the dispersion profiles 
of IScarpa et all \200H ) with a MOG dispersion profile and 
showed that the flattening of the profiles was consistent with 
Newtonian gravity, although the authors di d omit one data 
point from two of the clusters analysed. The IMoffat fc Tothl 



(2008) MOG model differs from Newtonian gravity for very 
massive objects (e.g. galaxy clusters and elliptical galax- 
ies) but becomes Newtonian for masses below a few times 
10 6 Mq , which is the mass range of GCs. The MOG fit proves 
to be an equivalently good fit to these data, indicating no 
deviation from Newtonian gravity in these clusters, and a 
Gaus sian with a fl a t tail has no physical significance. Be- 
cause IScarpa et all (|2007 | ) overlaid their data with a Gaus- 
sian with a flat tail, and IMoffat fc Tothl (|2008h overlaid the 
same data with a MOG model (which becomes Newtonian 
at low mass), and both are a similarly good fit, neither DM 
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r s = 2.3±1.2pc 
Central dispersion = 5.0±0.9kms 
Total mass = (9.0+"),10 4 M s 
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Figure 3. Velocity dispersion profiles of each cluster. The best fitting Plummcr] dl91 ll ) model is overplotted, and the fitted paramete 
as well as the calculated total cluster mass is shown. 



nor a modified version of gravity are required to explain 
their data, and our study confirms this. 



3.3 Systemic Velocities 

In addition to the velocity dispersions, the systemic veloc- 
ity of each cluster was measured (using the Monte Carlo 
method desc ribed above) and compared to similar data in 
iHarrid |l996t ). the most comprehensive catalogue of GC pa- 
rameters. This comparison is summarised in Table [2] Note 
especially M53, for which our cal c ulated velocity is signifi- 
cantly different to that of Irlarrisl l| 19961 ). with the discrep- 
ancy w e ll out side the error bars of both measurements. Since 
Harri s! (1199a) q uotes systemic velocity data for M53 from 
Webbinkl(|l98ll ). who examined a total of 12 stars, our sys- 



temic velocity calculation based on 180 stars is a more reli- 
able measurement for this cluster. 



3.4 Mass-to-Light Ratio 

One important indicator of large quantities of DM in a dy- 
namical system is a large mass-to-light ratio (i.e. M/Ly 2>1). 



Table 2. Comparisons between the radial velocities of each clus- 
ter in the lHarrij (|l996h catalogue and those from this survey. 



Cluster 


V, ("Harris 19961 


Vr- (this paper) 


M22 


-148.9±0.4 


-144.86±0.34 


M30 


-181.9±0.5 


-184.40±0.20 


M53 


-79.1±4.1 


-62.80±0.31 


M68 


-94.3±0.4 


-94.93±0.26 



Furthermore, this DM will cause the stars to have large ac- 
celerations, and therefore inherently higher maximal veloci- 
ties, and hence a larger velocity dispersion. An obvious way 
to determine whether or not a pressure-supported object like 
a GC is dark matter dominated is to measure its M/Ly. Of 
course, this require s a surface brightnes s profile, for which we 
have used those bv lTrager et al\ (| 19951 ). These profiles were 
converted to solar luminosities per square parsec. The den- 
sity p rofile was then calculated using (again, see iDeiongliel 
Il987l ): 



p(r) 



Mtot 



(rl+r 2 ) 2 
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The M/Ly profiles are shown in Figure [4] along with the 
M/Ly ratios. The vertical line shows the value of r s . Be- 
cause the central regions of GCs are highly concentrated, 
particularly in post-core-collapsed clusters, it is difficult to 
accurately determine the mass in the core. The M/Lv ratios 
hav e, therefore, only b e en ca lculated for radii greater than 
r s . lEvstigneeva et all |2007) showed that ultra-compact 
dwarf galaxies (UCDs) follow the same luminosity - veloc- 
ity dispersion relation as old GCs. Since it has also been 
shown that UCDs with a dynamical M/Ly up to 5 do not 
requi re dark matter to explain these mass-to-lig ht ratios 
(e.g. Ha segan et al. 1 120051 ; lEvstigneeva et all 120071 ). we see 
no requirement for any of our cluster s, apart from M53, to 
contain any dark matter component. lHasegan et all (|2005T ) 
determined that a M/Lv > 6 for UCDs may indicate some 
dark matter content, and since we estimate M/Ly ~ 6.7 for 
M53, this may indicate a small dark matter component in 
this cluster. Note that if we take the lower limit of our un- 
certainties for the mass-to-light ratio of M53 we reach the 
regime where dark matter is not required, and none of our 
clusters show M/Ly 3> 1, indicating DM is not dominant. 

Note that we have produced M/Ly profiles rather than 
simply deriving a single M/Ly value from the central ve- 
locity dispersion and central surface brightness. We sug- 
gest that determining M/Ly profiles is preferable because 
it describes the M/Ly of the entire cluster outside the core, 
rather than just its core. This is especially true for post-core- 
collapse clusters, such as M30, in which the core masses are 
highly uncertain. Furthermore, M30 is the o nly cluster in 
our s ample known to have a collapsed core (e.g. lTrager et all 
1995). This cluster also has a maximum in its M/Ly profile 
at about the scale radius r s , which may be due to mass segre- 
gation. The most massive stars are known to f all toward the 
core during the evolution of the cluster (see ISpitzerl 1 19851 . 
for a review of mass segregation processes in GCs). The 
somewhat less pronounced maximum in the M/Ly profile of 
M68 may be the first evidence that this cluster is currently 
undergoing core collapse. It has been argued that a central 
concentration parameter c m 2 may in dicate a collapsed cor e 
in certain GCs, and for M68 c = 1.91 jvan den Berghll 19951) . 
Note that, despite our calculation of the M/Ly ratio be- 
ing restricted to R > r s , the profiles should be accurate for 
R > r c , the core radius, or R > r s /2. Therefore, we argue 
that the peak in M/Ly is real, although a larger sample of 
cluster members for R < r 3 will be required to test this and, 
therefore, to show whether M68 is truly undergoing core 
collapse. 



4 CONCLUSIONS 

We have studied four GCs to determine their velocity dis- 
persion profiles, and found that neither DM nor modified 
gravity are required to reconcile the kinematic properties 
of these particular clusters. M53 may contain a small DM 
component, indicated by a M/Ly > 6, but it is not dom- 
inant. Within the uncertainties, the M/Ly ratio of M53 is 
still consistent with little, or no, DM content. The dynam- 
ics of all four c l usters are well described by a purely an- 
alytic iPlummerl l|l91ll ) model, which indicates that New- 
tonian gravity adequately describes their velocity disper- 
sions, and shows no breakdown of Newtonian gravity at 



ao ~ 1.2 x 10 'ms~ , as has been claimed in previous 
studies. 

Furthermore, the Plummer model was used to deter- 
mine the total mass, scale radius (the radius at which half 
the mass is contained), and the M/Ly profile for each clus- 
ter. We find that none of our clusters have M/Ly 2> 1, 
another indication that DM is not dominant. Within the 
uncertainties, our est imated cluster masses all match those 
in the literature fe.g. iMeziane fc Colin 1996 ), and the same 
is true for the M/Ly ratios (e.g. lPrvor fc Mevlanlll993h . We 
have produced M/Ly profiles, rather than quoting a single 
value based on the central velocity dispersion and central 
surface brightness. This method is preferred because it de- 
scribes the M/Ly of the entire cluster, rather than only its 
core; this is particularly important for post-core-collapsed 
GCs for which the central mass is highly uncertain. 

Because the only known of our clusters to have a col- 
lapsed core (M30) shows a peak in its M/Ly profile at 
approximately r s , this may be an indication of the re- 
sulting mass segregation. If this is the case, then we may 
have the first evidence that M68, which also exhibits a 
M/Ly peak near r s a nd has a central concentration of 1.91 
llvan den Bcrgh 1995), is currently undergoing core collapse. 
Further observations are required, both photometric and 
kinematic, to confirm this. 

Another important result from this study is the mea- 
sured rotations of our clusters. Of the four clusters studied 
here, one exhibits an obvious rotation, namely M22, with 
an axis of rotation approximately North-South. M68 may 
be rotating as well and additional data are required to con- 
firm this. M53 and M30 do not show any rotation to the 
level of 0.5kms _1 and 0.8kms _1 respectively. 

While these results are strongly indicative of the current 
picture of globular clusters being dark matter poor, and their 
dynamics can be explained by standard Newtonian theory, 
we are currently undertaking robust dynamical modelling of 
GC systems to fully address the question of the influence of 
non-Newtonian physics in explaining our observations. 
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Figure 4. Mass-to-light profiles for each cluster (thick curve). The thin curves show the uncertainties on the M/Lv, calculated via the 
1 a difference of r 3 from its calculated value. The vertical line indicates the value of r s . The M/Lv value given is only calculated for radii 
greater than r s because the mass estimates are highly uncertain in the cores of GCs, particularly if they are post-core-collapsed clusters. 
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